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I. INTRODUCTION 



Driven by the heavy-ion programs at GSI, CERN SPS, 
RHIC and LHC there is strong interest in the properties 
of strongly interacting matter at extreme temperatures 
and baryon densities. The access to the phase diagram 
of QCD is hampered by the fact that lattice simula- 
tions at finite density suffer from the sign problem 
However, in recent years impressive progress has been 
made in the understanding of the phase diagram of QCD 
within QCD effective models such as the Polyakov-loop 
extended Nambu-Jona-Lasinio (PNJL), see e.g. [H, Hf, 
and the Polyakov-loop extended quark-meson (PQM) 
model, see [6|. In these models the information about 
the confining glue sector of the theory is incorporated in 
form of a Polyakov-loop potential that is extracted from 
pure Yang-Mills lattice simulations. The matter sector of 
these models has been studied in detail also beyond the 
mean-field level by taking into account the quark-meson 
quantum fluctuations mostly within a functional renor- 
malization group (FRG) approach, for recent studies see 

an. 

However, the most difficult problem within these 
models is the question of how to embed the quantum 
back-reaction of the matter sector to the gluonic sec- 
tor. This problem has effectively been resolved in Q 
where the change of Aqcd, and hence the confinement- 
deconfinement transition temperature To, in the pres- 
ence of dynamical quarks has been computed within hard 
thermal and/or hard dense loop (HTL/HDL) perturba- 
tion theory. This leads to a flavor and chemical poten- 
tial dependence of the transition temperature To, and 
is a qualitatively viable procedure: the change of the 
non-perturbative scale Aqcd can be very well estimated 
within perturbation theory as has been shown and con- 
firmed in many computations at zero-temperature and 
finite-temperature QCD. Moreover, this perturbative es- 
timate in [f| has been confirmed by first principle QCD 
computations with the functional renormalizati on g roup 
at real and imaginary chemical potential in ha. On 



the other hand, very recently the /i-dependencc of To 
has been estimated by constraining PNJL results with 
those in the statistical model ff|. This links the model 



parameter To(iV/,/x) under certain weak assumptions to 
experimental data about the chemical freeze-out curve 
and nicely confirms the theoretical prediction in Q . 

This shows that even though being QCD effective mod- 
els, they provide valuable information about the phase 
diagram, and in particular help to exclude certain scenar- 
ios. Moreover, the models can be understood as specific 
controlled approximations to full dynamical QCD. Most 
directly this is realized within the above-mentioned FRG 
approach put forward in 0, [hJ If2l - fl4| . This link enables 
us to systematically extend the PNJL/PQM models to- 
wards full dynamical QCD. 

In the present work we take an important step towards 
full dynamical QCD and present the first computation 
with fully dynamical matter sector in the PQM model at 
finite density. An interesting first FRG computation at 
vanishing density has been put forward in [15|. There, 
however, the back-coupling to the glue sector has been 
neglected. Here, we include the back-reaction of the mat- 
ter sector as described above. The uncertainty can be 
estimated by comparing it to the QCD computation in 
0, LT3] ■ Within the present PQM approach beyond mean 
field we provide results on the phase boundaries for the 
confinement-deconfinement and the chiral transitions as 
well as on thermodynamic quantities such as the pressure, 
quark number density and quark number susceptibility. 

The outline of the paper is as follows: in the next 
section we discuss the back-reaction of the matter fluc- 
tuations to the QCD Yang-Mills sector in more detail. 
In Sec. IIII1 the Polyakov-loop and its potential is intro- 
duced and coupled to the quark-meson model which then 
defines the Polyakov-quark-meson model. At first, the 
grand potential of this model is evaluated in mean-field 
approximation. Afterwards, the flow equation for the 
grand potential of this model is derived in Sec. IIVI where 
also the choice of model parameters is discussed. Sec. Mis 
devoted to the phase structure and some thermodynam- 
ical applications, in particular we evaluate the pressure, 
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quark number density and quark number susceptibility. 
The phase structure, i.e., the chiral and confinement- 
deconfinement phase transition is explored in detail. Sub- 
sequently, the influence of the Polyakov loop on the ther- 
modynamics is investigated and in Sec. I VII concluding 
remarks are drawn. 



II. THE PQM MODEL BEYOND MEAN FIELD 

The Polyakov-quark-meson (PQM) model was stud- 
ied in a mean-field approximation for the matter sector 
in . In this model, the coupling between the pure glue 
sector, realized by an effective Polyakov- loop potential 
U, and the quark- meson matter sector is provided by the 
fermionic determinant. However, in Q an important step 
beyond the mean-field approximation was already intro- 
duced. More specifically, the pure glue potential U was 
adjusted on the basis of phenomenological arguments: 
the phase transition temperature To in the potential U 
relates to the dynamical scale in QCD, Aqcd- If the 
dynamics of the quarks is switched on, this scale is low- 
ered. In we have provided a flavor and quark chemi- 
cal potential dependent transition temperature Tq(Nj, /i) 
which was estimated from hard thermal and hard dense 
loop considerations. Here, we shall elaborate this argu- 
ment in more detail, also in order to show that no double 
counting is involved in the computation and to tighten 
the link to first principle QCD. We also stress that the 
above phenomenological argument works well at large 
chemical potential fi as well as at small chemical poten- 
tial. For vanishing chemical potential we also compare 
our results with a recent dynamical QCD calculation Q . 
This allows us to qualitatively adjust the phase transition 
lines at small chemical potential and small temperature, 
and gives us access to information about a possible crit- 
ical endpoint as well as the size of a possible quarkyonic 
matter region at large chemical potential in the QCD 
phase diagram as suggested in [161 ] . 

What is then left is to include the quantum and ther- 
mal fluctuations of the matter sector in the presence of 
the non-trivial Polyakov loop. This is achieved with func- 
tional renormalization group (FRG) methods and is done 
in Sec. IIV1 for reviews see e.g. [l?], [Hi- I n combination, 
this allows for a Polyakov-loop effective model computa- 
tion which already takes into account the full dynamics 
of QCD phenomenologically. 

We proceed with the explanation of the Nf and /i de- 
pendence of the transition temperature To. In a full 
two-flavor QCD computation was put forward for vanish- 
ing and imaginary chemical potential. The related func- 
tional RG equation for the effective action is provided 
in a simple diagrammatic form in Fig. [TJ The first two 
loops in Fig. [TJ stand for quantum fluctuations (gluons 
and ghosts) in the pure glue sector that, e.g., generate 
the Polyakov-loop potential. The third loop represents 
the quark fluctuations, and the last term encodes mcsonic 
fluctuations generated by dynamical hadronization (l8l - 



FIG. 1: Functional QCD flow for the effective action: The 
lines denote the corresponding full field dependent propaga- 
tors. Crosses denote the cut-off insertion dtR. 

[20j . The crosses denote the cut-off insertion which re- 
stricts the loop momenta to that about the cut-off scale, 
and (j> stands for all fields. For more details see [H, [1(3, EH . 

We would like to emphasize that the above equation 
for the effective action is fully coupled. An important ex- 
ample is the gluonic (and indirectly the ghost) propaga- 
tion that is modified in the presence of dynamical quarks. 
This aspect is visualized in Fig. [2] where the quark contri- 




FIG. 2: Quark contributions to the flow of the gluon propa- 
gator 

bution to the gluonic propagator flow is shown. In turn, 
there are gluonic contributions to quark and meson cor- 
relation functions, and the system is highly non-linear. 

Another important aspect is that the pure Yang-Mills 
(YM) Polyakov-loop potential can be derived from the 
pure Yang-Mills flow as depicted in Fig. [3] This simply 



d t T k [(j)] = 




FIG. 3: The pure Yang-Mills flow for the effective action. 

amounts to dropping the matter loops in Fig. JTI The 
pure Yang-Mills computation has been done in [21|, and 
the related To agrees quantitatively with lattice predic- 
tions which are used to adjust the To parameter in the 
Polyakov-loop potential U. These findings also elucidate 
that To is nothing else but the dynamical YM-scale Aym 
inherent in the ghost and gluon propagators. This sup- 
ports the quantitative accuracy of such an approach. 

In turn, the quantum dynamics of the quark-meson 
model is obtained by dropping the Yang-Mills diagrams 
in Fig. [TJ, that is the gluon and ghost loop. This has been 
studied extensively in the literature, for recent works 
see e.g. 0, H|- In summary, it is this simple additive 
structure for the different contributions that allows to 
systematically improve quark-meson models towards full 
QCD within an FRG approach. Moreover, it also al- 
lows us to directly use full QCD information within 
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the PQM/PNJL models. This reduces the parameter- 
dependence of these models, and qualitatively enhances 
the predictive power. 

An important example for this structure is the 
Polyakov-loop potential in these models: the Polyakov- 
loop potential U in full dynamical QCD still comes from 
the Yang-Mills contributions in Fig. [T] The related 
Polyakov-loop potential in QCD can be computed from 
the full dynamical QCD computation in (9|, and receives 
a flavor and chemical potential dependence via Fig. [2] It 
is clear that the quark contributions to the gluon propa- 
gation (and ghost propagation) shown in Fig. flower the 
dynamical scale in the gluon and ghost propagation, and 
hence To . In turn, with constant To without a flavor and 
/^-dependence, this contribution is estimated by the pure 
Yang-Mills results. We conclude that the hard thermal 
and hard dense loop adjustment of To is necessary for 
the full inclusion of quark effects to the Polyakov-loop 
potential. In particular, the above contributions are not 
included in the fermionic determinant and thus no double 
counting is involved. 

In summary, the inclusion of the Nf and ^-dependence 
in T) as well as the dynamics of the quark-meson sector 
of the model in the presence of a non-trivial Polyakov- 
loop expectation value gives us a good qualitative control 
about the full QCD dynamics, and the link is tightened 
by comparing T (N f , fx) with QCD results in @,[h|. The 
direct use of the full QCD Polyakov-loop potential will 
be discussed elsewhere. 



The isoscalar-scalar a field and the three isovector- 
pscudoscalar pion fields 7? together form a chiral vector 
field 4>. Without the explicit symmetry breaking term 
c in the mesonic potential the Lagrangian is invariant 
under global chiral SU(2)l x SU{2)r rotations. 

The parameters of the model defined in Eq. (UJ arc 
fixed at the low energy physics at vanishing temperature. 



A. Polyakov-loop potential 

In recent years, several Polyakov-loop potentials have 
been proposed [22|, HH , for comparison see [24| . Their 
functional form is motivated by the underlying QCD 
symmetries in the pure gauge limit and they all repro- 
duce a first-order transition at T c ~ 270 MeV for N c = 3 
colors in this limit. 

In this work we use a polynomial ansatz for the 
Polyakov-loop potential 



Wpoly _ b 2 (T) 



' J 3 / S 3 



$3 + $ 



3 ^ ' 4 '$$) 2 , (4) 



T 4 2 6 v ' ' ' 4 

with the temperature-dependent coefficient 

b 2 (T)=a + a 1 (^ya 2 ^ 2 + a,^ 3 . (5) 

The parameters of Eqs. @ and §5§ are fitted to the 
pure gauge lattice data with 



III. THE POLYAKOV-QUARK-MESON MODEL 

The Polyakov-quark-meson model as put forward in 
already provides a good approximation of full QCD at 
low energies and temperatures and not too high densities. 
Its classical Euclidean action with Nf = 2 light quarks 
q = (u, d) and N c = 3 color degrees of freedom reads 

d 4 xl q(p + h(a + z 75 T7?)) q + \{d^) 2 



+ ~(d^) 2 + U(a, tF)+W(M) 



(1) 



where fl(&) = 7^9^ — ig^foA (^) with gauge coupling 
g and h is the Yukawa coupling between mesons and 
quarks. The temporal component of the gauge field is 
linked to the order parameter of quark confinement (in 
pure Yang-Mills), the Polyakov-loop variable $, 

$(x) = -^/trPexp^ 5 ^dr^o(x,r)^\ , (2) 

where the color trace tr is in the fundamental representa- 
tion, V denotes path ordering, and /3 = 1/T, the inverse 
temperature. The purely mesonic potential is defined as 

U(a,7r) = ^(a 2 +tt 2 -v 2 ) 2 -ca . (3) 



a = 6.75 , ai = -1.95 , a 2 = 2.625 , a 3 = -7.44 (6) 
and 



0.75 



7.5 



(7) 



The above potential reproduces well the equation of state 
and the Polyakov-loop expectation value, in particular 
around the transition [24j. The parameter To = 270 MeV 
corresponds to the transition temperature in the pure 
YM theory. 

As motivated in the previous section, it is left to fix To 
in the presence of dynamical quarks. Since it is directly 
linked to the dynamical mass-scale Aqcd this parameter 
necessarily has a flavor and chemical potential depen- 
dence in full dynamical QCD and To — > To(Nf,fi). 

In the present work we use perturbativc relations for 
fixing the relative scales (fj HH . The results compare well 
to the full Nf = 2 QCD computation in the chiral limit 
in 0, . The latter thus allows for an error estimate of 
the present procedure. The one-loop /3-function of QCD 
with massless quarks is given by 



/3(a) = -ba 2 , 



with the coefficient 
b(N f ) 



6tt 



(lliV c - 27V/) . 



(8) 



(9) 
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Here, we have assumed a RG scheme that minimizes (part 
of) the higher-order effects. At leading order the corre- 
sponding gauge coupling is given by 



a{p) = 



1 + a b{N f ) ln(p/A) 



+ O(og), 



(10) 



with ao = a(A) at some UV-scale A. The scale 
Aqcd = Aexp(— l/(ao6)) corresponds to the Landau 
pole of Eq. (fTOj), 

The temperature dependence of the coupling is also 
governed by Eq. (fill with the identification p ~ T . This 
yields the relation g 



T (N f ) =fe- 1/{a ° b(N f )) , 



(11) 



where T and ao are free parameters. Eq. (1111) allows 
us to determine the TVy-dependence of the critical tem- 
perature T (Nf). Analogously to @ we choose T to be 
the r-scale, T = T T = 1.77 GeV. This constitutes a rea- 
sonable UV scale for the mean-field model. Then the 
pure Yang-Mills input, T (N f = 0) = 270 MeV, leads 
to ao = 0.304. In the present work we shall stick to 
these values. In addition to the arguments given in Q, 
the ratio Tq/T x in the chiral limit compares well with 
that computed in the full two-flavor QCD calculation in 
Q . Table U summarizes the A/-dependent critical tem- 
perature To in the Polyakov-loop potential for massless 
flavors: 







2 + 1 



To [MeV] || 270 | 240 | 208 | 187 | 178 

TABLE I: Critical Polyakov-loop temperature To for Nf mass- 
less flavors. 



Massive flavors lead to suppression factors of the order 
Tq/(Tq + to 2 ) in the /?- function. For 2 + 1 flavors and 
a current strange quark mass m s ~ 150 MeV we obtain 
To (2 + 1) = 187 MeV. We estimate the systematic error 
for T (N f ) being of the order j^o MeV related to the 
scale matching of the present PQM computation with 
the QCD computation in the chiral limit in Note, 
however, that the link to QCD qualitatively improves 
the error estimate in comparison to the estimate done in 

!• 

As argued in the last section, in addition to the fla- 
vor dependence of To we introduce a chemical poten- 
tial dependence via a ^-dependent running coupling fo, 
which should push the confinemcnt-dcconfincmcnt tran- 
sition temperature down close to the chiral transition 
line. This can be achieved by defining 



with 



T (N f ,fi) = T T e- 1/iaob{N f'^ 



b{N f ,ri=b(N f ) 



'(7 TV) 



(12) 



(13) 



The factor 7 is a parameter governing the curvature of 
T ((i) and 6 M ~ ^§-Nf as in 0. As for the iV^-dependence 
the ^-dependence in Eq. (fl"2| compares well to that found 
m QCD Based on the results there we estimate 

the systematic error with 0.7 < 7 < 1, and we shall 
investigate the 7-dependence of our results in Sec. [Vj 



B. Grand Potential in Mean-Field Approximation 

All thermodynamic properties of the PQM model fol- 
low from the grand potential. It is a function of the 
temperature and one quark chemical potential since we 
consider the SU (2) /-symmetric case in this work and set 
/i = fj, u = [id, 

In the mean-field approximation certain quantum and 
thermal fluctuations in the path integral representation 
of the grand potential are neglected. The mesonic quan- 
tum fields are replaced by their corresponding classi- 
cal expectation values and only the integration over the 
quark loop is performed which is modified by constant 
gluon background fields in the PQM model [2f|. The 
final potential in mean-field approximation reads 

fl MF = n qq (a, $, $) + U(a, 0) + $) (14) 

and consists of the quark contribution including the 
Polyakov-loop variables 



fi ? , = -2N f 



Xg-^p-^/T + e -3(E p -u)/T 



In 



1 + 3(# + $ e -( £ p+^)/ T )e-( iJ p+^/ T 



+ e -3(E p +n)/T I 



(15) 



with the quark/antiquark singlc-quasiparticlc energies 



- . ^2 



lp" + to 2 , and the constituent quark mass m q = 

ha. The purely mesonic potential U is given by Eq. ([3]) 
and the effective Polyakov-loop potential U, e.g., by 
Eq. ([4]). Details of the potential derivation can be found 
in [6|. The quark contribution involves a divergent vac- 
uum term which can be regularized. As shown in (27l . |28| 
this term is important and modifies the underlying ther- 
modynamics. Since this term upgrades the standard 
mean-field approximation it is neglected here whereas it 
is included in the full RG approach. 

The solution of the corresponding equations of motion 
are obtained by minimizing the thermodynamic potential 
with respect to the three mean fields a, $ and i.e., 



on 



MF 



MF 



on 



MF 



da 



9$ 



0. 



(16) 

The solutions to Eq. (jTB")) provide the chiral (a) and 
Polyakov-loop expectation values ($} and ($) as func- 
tions of the temperature and quark chemical potential. 
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IV. FLOW EQUATION 

The non-perturbative FRG method has a wide range 
of applicability. In the context of equilibrium statistical 
physics it represents a very efficient way to describe crit- 
ical phenomena and in particular phase transitions. One 
particular formulation of RG flows is based on the con- 
cept of the effective average action Tk where k denotes 
a RG momentum scale, for reviews see e.g. [ItI EH- For 
a system with bosonic (ip) and fermionic fields (ip), the 
variation of with the RG scale (t = In k) is governed 
by the flow equation 

d t T k [ipM = ^Tr jdA^r^WH^Mj) 1 } 

- Tr id t R k>F (if' 2) [ip, V] + Rh,f) ~* | 

(17) 

where 1^ denotes the ith (jth) derivative of Yk with 
respect to the <p (ip) fields. The trace involves a d- 
dimensional momentum integration and a summation 
over all inner spaces (flavor, color and/or Dirac). The 
regulator functions R^i and their derivatives implement 
Wilson's idea of integrating successively over narrow mo- 
mentum shells and ensure that the flow equation is both 
infrared and ultraviolet finite. We employ the optimized 
regulator functions [29| which depend only on the spatial 
components of the momentum and for bosonic fields is 
given by 

i? fe , B = (fc 2 -/)e(i-g) (is) 

whereas the fermionic regulator for the PQM model reads 

r*s=*(J%- -g) • (19) 

I 



Further details and the finite temperature and density 
generalization of the flow equation can be found, e.g., in 
Refs. jllUEi- 

In the PQM model with the classical action, Eq. (JTJ, 
the flow Eq. (|17p for the quark-meson sector also de- 
pends on the Polyakov loop via the Dirac term. Here we 
will take into account the full dependence of the quark 
and meson fluctuations on a general constant $, $ back- 
ground. In summary, this corresponds to the truncation 



r fc = / d 4 x {i; (p + f i lQ + ih(a + n 5 TTr)) ip (20) 



where Clk[c, 7?, $] is the full effective or grand poten- 
tial in a general background $ and $. It is the quan- 
tum analogue of the mean- field potential J^mf, given in 
Eq. (p~4|) and, in the present truncation, carries the full 
fc-dependence of the effective action T^. The approxima- 
tion Eq. (|20p has also been used for the matter sector 
in @ and in the first PQM study beyond mean field at 
vanishing density [l5j . Finally, we are interested in flk=o 
evaluated at the solutions of the quantum equations of 
motion. This corresponds to the thermodynamic poten- 
tial where from all temperature, chemical potential and 
field derivatives follow. The spatial gauge field is set to 
zero while we keep the temporal component as a con- 
stant mean-field background. With this truncation the 
flow equation for the effective potential in leading order 
derivative expansion reads 



9 t fi fc [a,7?,$,$] = 



k* 



12tt 2 



-^-coth ( — - 
E n \2T 



— coth [ — - 
E„ \2T 



2ug 

E n 



{1 - iV 9 (T, M ; $, !>) - 7V g -(7>; $, $)} 



(21) 



with the Polyakov- loop enhanced quark/antiquark occupation numbers 

1 + 2$>e {E «-^/ T + $ e 2(E,-M)/T 



N q (T, 



I 



and N q (T, ,i; $, $) = N g (T, -fi; $, $) , (22) 



see also the QCD study [9J at vanishin g an d imaginary 
chemical potential, and the PQM study [15| at vanishing 
chemical potential. In the present work we explore the 
full phase diagram at real chemical potential. The num- 
ber of internal quark degrees of freedom is denoted by 
v q = 2N c Nf = 12. This equation describes the flow of 



the full quark and mesonic subsystem modified by a con- 
stant background field. The grand potential also depends 
on the Polyakov-loop variables and the expectation value 
of the square of the chiral 4-componcnt field <fi 2 = a 2 + 7? 2 
which coincides with (a) 2 since {7?} 2 = 0. 

The quasi-particle energies for i = q,<J,-K are given by 
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100 150 200 250 300 100 150 200 250 300 

T [MeV] T [MeV] 

FIG. 4: The normalized chiral condensate a/ao and the Polyakov-loop variable $ (left panel) and the corresponding temperature 
derivatives (right panel) as a function of the temperature for vanishing chemical potential. The Polyakov-loop parameter 
To = 208 MeV encodes the back-reaction to the glue sector of the quark dynamics in the presence of two flavors. 



Ei = \Jk 2 + mf with the corresponding squared con- 
stituent quark and meson masses respectively 

m\ = h 2 <? , ml = 2n' k + , ml = 2n' k . (23) 

Primes denote the 2 -derivative of the grand potential, 
i.e., Q' k :=dn k /d<j) 2 . 

In the limit of vanishing background fields, i.e., when 
<£>, $ — > 1, the extended occupation numbers simplify to 
the usual Fcrmi-Dirac distribution functions for quarks 
and antiquarks 



N q (T,fi;l,l) = y 



1 



(24) 



exp((£, - M )/T) 

A^T^M) = — 1 (25) 

l + exp((£'g +fJ,)/T) 

and the flow of the quark-meson model is recovered [H, 

The flow equation (|2T|) constitutes a set of coupled, 
highly non-linear partial differential equations that can- 
not be solved with analytical methods since the right 
hand side of (j2l~j) depends on derivatives of the unknown 
potential tt k . For the sake of full quantitative precision 
and in order to facilitate the access to the potential first- 
order phase transition region in the QCD phase diagram 
we do not resort to Taylor expansions but rather com- 
pute the full numerical solution of the effective potential 
Ofc[cr, 7r, $, $]. This extends the analyses of the full effec- 
tive potential in the quark- meson model, see to the 
present case. The solution of Eq. pip yields the ther- 
modynamic potential fi G ff [er, 7?, $, <I>] = Q k=0 {<j, ir, $, <I>] 
for the Polyakov-extended quark-meson model in the in- 
frared. The expectation values of the fields are deter- 
mined by the quantum equations of motion, 



dn cff dn cS d£l cS 



da 



d<S> 



<9$ 



0, (26) 



in analogy to the mean- field analysis in Eq. (|16l) . 



A full study including algorithmic differentiation (AD) 
techniques [36[ will be presented elsewhere [13] ■ In the 
present work we shall simply evaluate the effective po- 
tential 0^ on the solutions $(ct),$((t) of the mean field 
EoMs. We shall argue that this is already a quantita- 
tively reliable approximation to the full solution: note 
that the present truncation introduces a cut-off, and 
hence implicitly a momentum-dependent fermion propa- 
gator. For large temperatures T in the deconfined regime 
the fermion propagator is cut-off by the Matsubara mass 
7rT. At small temperatures the theory is chirally bro- 
ken and the fermion propagator exhibits a mass of the 
order of the chiral scale ~ 360 MeV. Moreover, at suffi- 
ciently large densities the fermionic fluctuations are al- 
ready well-described by the one loop determinant. In 
summary, the fluctuation dependence of the Polyakov 
loop can be treated as a perturbation in the whole phase 
diagram. Accordingly, the back-reaction of the Polyakov 
loop beyond mean field to the matter fluctuations is sup- 
pressed for all T and [i. Indeed, in comparison to the 
full solution the minimum for a varies within ±3 MeV 
whereas <E> are naturally more sensitive to this approx- 
imation and vary within ±20 MeV. Note, that the above 
structure provides further non-trivial reliability for the 
Polyakov-loop extended models. Seen as an expansion 
towards full QCD they partially allow for perturbativc 
arguments about the mean-field analysis. The full com- 
putation and a more detailed analysis of this structure 
will be provided in [37j . The above structure also empha- 
sizes the crucial input: the back-reaction of the matter 
sector to the Polyakov loop via T (Nf,[j,). 

It remains to determine the initial effective action Ta, 
or, more precisely, the effective potential fl\ at the initial 
scale A = 950 MeV in the UV. First of all, consistency 
with the flow Eq. (f2~Tj) leads to a (relevant) term origi- 
nated in the fermionic loop, ^a°[ (J ' ™, ^ see Eq. (p?5j) . 
This term is relevant for the correct thermodynamics and 
also includes fermionic vacuum fluctuations. The field- 
dependent part of Q\ consists of a sum of the quark- 
meson potential Eq. ([3]) and the external glue input, the 
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Yang-Mills Polyakov-loop potential U. This yields 

ft A [o-,7f,$,$] = C/(cr,7r)+W($,$)+^[cr,7f,$,$] . (27) 

It is left to determine the model parameters in the 
quark-meson sector. As in the mean-field analysis, they 
are fixed to reproduce the physical quantities in the vac- 
uum such as the pion decay constant, /„• = 93 MeV, pion 
and sigma meson masses, = 138 MeV, m a = 500 
MeV, and the constituent quark mass m q = 298 MeV 
in the infrared. The explicit symmetry breaking term 
c = Jti^/tt is an external field and is scale independent. 



V. PHASE STRUCTURE 

We first discuss the order parameters at vanishing 
chemical potential and compare our findings with lattice 
results [H, HI] as well as continuum QCD results 0] . We 
evaluate the validity of our approach, in particular, the 
need for the inclusion of the back-reaction of the mat- 
ter sector via a Nf and //-dependence of the dynamical 
transition parameter To in the Polyakov-loop potential. 

Fig. 2] summarizes our findings for fi = and the dy- 
namical transition parameter To = 208 MeV. It shows 
in the left panel the chiral condensate, normalized with 
the vacuum value and the (degenerated) Polyakov-loop 
variable as a function of the temperature. In the right 
panel, the corresponding temperature derivatives of the 
order parameters are displayed. As one can see from 
Fig. 21 both widths in the temperature derivatives are 
comparable. For the chiral transition a slightly larger 
critical temperature, T c ~ 190 MeV as for the deconfine- 
ment transition, T c ~ 175 MeV, is found. This relates 
to the standard scenario that chiral symmetry restora- 
tion occurs at higher temperatures than deconfincment. 
We remark that the arguments for the standard scenario 
are strictly valid only for phase transitions rather than 
for smooth crossovers as in the present case. Moreover, 
the chiral and confincment-deconfincmcnt crossover tem- 
peratures agree within the respective widths. The above 
findings compare well with a two flavor QCD computa- 
tion in the chiral limit and also with recent lattice 
results [38|, |39[ . For the comparison with the latter one 
has to bear in mind that the definitions of the chiral order 
parameter are ambiguous which influences the crossover 
temperature. 

If we switch off the back-reaction of the matter sec- 
tor and stick to the pure Yang-Mills value of the tran- 
sition parameter, To = 270 MeV the situation changes 
quantitatively as can be seen from Fig. [5] the chiral 
and Polyakov derivatives peak at similar critical tem- 
peratures about T c ~ 220 MeV but the width of the 
Polyakov derivative is very broad suggesting a peak sub- 
structure. This broad structure does not compare well 
with the lattice predictions nor does the high crossover 
temperatures. This demonstrates the importance of the 
dynamics of To already at vanishing chemical potential. 



Now we extend our analysis to finite chemical poten- 
tial. Our findings are summarized in Fig. [5] The left 
panel shows, for comparison, the phase diagram with- 
out back-reaction of the matter sector to the glue sec- 
tor. The right panel shows the phase diagram with the 
full dynamics: the quark-meson dynamics in the pres- 
ence of the Polyakov-loop background is included with 
the flow Eq. ([2~Tj) whereas the back-reaction of the mat- 
ter sector to the glue sector is encoded in To(Nf,/i) as 
given in Eq. (fT2"j) . In the right panel of Fig. [5] we have 
used 7 = 0.85 and T (2,0) = 208 MeV which compares 
well to the QCD results in [To| for small densities. The 
shaded band corresponds to the width of d&/dT at 80% 
of its peak height. We observe, that the width of this 
band shrinks with increasing \x. Thus, the deconfinc- 
ment transition gets sharper at higher fi. Furthermore, 
we can unambiguously distinguish the peak positions of 
the chiral and deconfinement transition in the temper- 
ature derivatives. Even when we vary the parameter 7 
almost no differences emerge. Only for value of 7 close 
to one an ambiguity in the peak positions of the temper- 
ature derivatives arises, similar to the constant To-case 
discussed below. 

We see that the deconfinement transition line stays 
close to the chiral phase boundary. It has been specu- 
lated that the quarkyonic phase [16| is signaled by a re- 
gion with confinement and chiral symmetry Mean field 
computations in the PNJL and PQM model with con- 
stant To have supported this scenario. As has been shown 
in Q , this does not hold true if the dynamics of the tran- 
sition parameter To is taken into account. Here we have 
confirmed that in the fully dynamical PQM model the 
prediction in Q holds. 

The chiral first-order transition line arising at small 
temperatures and large chemical potentials terminates in 
a CEP which is located at fi C EP ~ 292 MeV, T CE p ~ 23 
MeV. As discussed in [|(| the precise location of the CEP 
depends on the chosen parameters in the vacuum, in par- 
ticular on the sigma meson mass. The back-bending of 
the first-order transition line in the phase diagram to- 
wards smaller chemical potential is typical for a FRG 
calculation Q. But very close to the //-axis the slope 
of the first-order line tends back to infinity similar to a 
mean- field treatment (not shown in the figure). This be- 
havior is also in agreement with the Clausius-Clapeyron 
relation according to which the transition line should hit 
the //-axis perpendicular. 

If we switch off the back-reaction of the matter fluc- 
tuations to the Yang-Mills sector and choose a constant 
To = 208 MeV, the above picture changes drastically. For 
finite chemical potential the Polyakov loops $ and $ start 
to deviate and the widths of their temperature derivatives 
increase over the whole phase diagram. The resulting 
phase diagram for a constant To is summarized in the left 
panel of Fig. [6l In the vicinity of the intersection point 
of the chiral transition and the deconfinement transition 
which are both smooth crossovers around /1 ^ 180 MeV a 
double peak structure in the corresponding temperature 
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FIG. 5: The normalized chiral condensate a/ao and the Polyakov-loop variable <E> (left panel) and the corresponding temperature 
derivatives (right panel) as a function of the temperature for vanishing chemical potential. The Polyakov-loop parameter To is 
set to its pure Yang-Mills value T = 270 MeV. 



derivatives of the Polyakov-loop variables occur. This 
hampers a unique identification of the transition point. 
In order to clarify this behavior we plot the maximum 
of the peak location of the T-derivatives in the phase 
diagram together with a band around 80% of the maxi- 
mum value. For larger chemical potential the peak loca- 
tions deviate strongly and a coincidence of the chiral and 
deconfinement transitions can be excluded. This brings 
back the chirally symmetric and confined region which 
has been connected to the quarkyonic phase. In sum- 
mary, we have shown that this signature is very sensitive 
to the correct implementation of the back-reaction of the 
matter sector to the glue sector, and is most likely not 
present. 

At high chemical potential and small temperatures, a 
first-order chiral phase transition takes place which ends 
in a critical endpoint of second order located at /xcep = 
293 MeV, T CEP ~ 32 MeV. 



A. Thermodynamics 

In general, the initial action given at the ultraviolet 
cutoff is independent of the temperature and chemical 
potential which limits the reliable calculation of thermo- 
dynamic quantities to a certain temperature and den- 
sity region. Due to the ultraviolet cutoff high thermal 
and quantum fluctuations are suppressed. As a conse- 
quence, cutoff-independent predictions can be obtained 
for temperatures T < A/8 which in our case yields an 
upper temperature bound of T ~ 120 MeV. In order to 
cure this cutoff remnant at high temperature one has to 
combine the FRG with the expected perturbative results. 
The inclusion of the missing high-momentum modes can 
be achieved in an effective way by adding to the orig- 
inal flow, Eq. (|2ip. a flow equation for an interacting 
Polyakov-loop quark system for scales k > A. Note, that 
an explicit gluon contribution to the flow equation is ne- 
glected here because the effective Polyakov-loop potential 
is fitted to reproduce the Stefan Boltzmann (SB) limit at 



high temperatures. In this way we employ the equation 
-N c N f k 4 



d k Q k A (T,fj,) = 



3ir 2 E n 



(28) 



[1 - A 9 ($, #; T, M ) - A g -($, $; 



with E a 



as previously defined. This equa- 



tion is integrated from k = oo to k = A and yields the 
UV-contribution f2^° (T,fi), which is then added to the 
grand potential resulting from the solution of the RG 
flow equation (|21j) . The divergent zero mode contribu- 
tion in Eq. (|28]l is neglected here. However, for vanishing 
quark masses this term represents an unobservable shift 
in the grand potential. 

In Fig. [7] the thermodynamic pressure normalized to 
the Stcfan-Boltzmann pressure for three different quark 
chemical potentials is shown MS irl function of tempera- 
ture. In the left panel of this figure, a fixed T = 208 MeV 
has been used while in the right panel the ^-corrections 
are taken into account. The inset displays the pressure 
for /i = 290 MeV which is close to the critical endpoint in 
the phase diagram. For small temperatures the pressure 
decreases and has a kink at the critical temperature due 
to the first-order transition. Without the back-reaction 
of the matter fluctuations to the Yang-Mills sector a simi- 
lar behavior in all three curves is observed. In the vicinity 
of the chiral transition the pressure increases due to the 
melting of the quark masses and saturates at about 80% 
of the corresponding ideal gas limit which reads for Nf 
massless quarks and (N% — 1) massless gluons 



Psb 

rp4 



N f N c 
6 



7tt2 
30 



I 

2^2 



+(iV c 2 -l) 



45 ' 

(29) 

Including the back-reaction of the matter sector via the 
inclusion of Tq((i) changes the thermodynamics at larger 
chemical potential. The pressure increases much faster 
and saturates at /x = 290 MeV at 95% of the Stcfan- 
Boltzmann limit. 
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FIG. 6: Chiral and deconfinement phase diagram for a constant To = 208 MeV (left panel) and for To(/t) with 7 = 0.85 (right 
panel). The (grey) band corresponds to the width of dQ/dT at 80% of its peak height. Close to the intersection point of 
the chiral transition and the deconfinement transition at mid chemical potential a double peak structure in the temperature 
derivative of the Polyakov-loop variables emerges. The (green) dashed line in this region follows the highest peak. 




FIG. 7: Pressure normalized to the Stefan-Boltzmann pressure for a constant To — 208 MeV (left panel) and with /^-corrections 
(right panel) for three different chemical potentials. The CEP is located approximately at /1 = 293 MeV. The insets show the 
pressure at fj, = 290 MeV for small temperatures. 



A similar trend is seen in the entropy, Fig. El and 
quark number density, Fig. |H1 if the //-corrections are 
taken into account. The entropy density decreases for 
small temperatures at fi = 290 MeV since the number 
of active degrees of freedom decreases when approach- 
ing the first-order transition from below. At the transi- 
tion the entropy jumps. The bump around T ~ 90 MeV 
(left panel) is a remnant from the smooth chiral crossover 
transition. This effect is completely washed out when the 
//-corrections are included (right panel). Similar to the 
findings for the pressure these corrections become more 
significant at larger chemical potential. 

This also appears in the quark number density n q = 
—dCl/dfj, which is plotted in Fig. [9] For comparison the 
corresponding SB-limits (dashed lines) are also shown 
in this figure. The quark density approaches the SB- 
limit always from below. Without the /i-corrections the 
Polyakov loop suppresses the quark densities for chemi- 
cal potential larger than the intersection point of the chi- 
ral and deconfinement transition in the phase diagram 
Fig. [5] With the Tq(/j,) corrections both transitions coin- 



cide over the whole phase diagram and as a consequence 
the quark number density approaches much faster the 
SB-limit (right panel of Fig. [5). 



In Fig.[T0]the scaled quark number density (left panel) 
and the corresponding scaled quark number susceptibility 
(right panel) for three different temperature slices around 
the critical endpoint (T CEP , T CEP ±5 MeV) as a function 
of the quark chemical potential are collected. In this fig- 
ure the //-corrections in To are omitted while in Fig. [TT] 
they are taken into account. Due to the chiral critical 
endpoint which is a second-order transition the suscepti- 
bility diverges with a certain power law [4l| . There are 
no strong modifications in the structure of the suscepti- 
bility divergence if the back-reaction of the matter sector 
is taken into account or not. As a consequence it seems 
that the size of the critical region around the CEP is not 
strongly modified by these fluctuations. The only differ- 
ence is that including the /t-corrections the peak height of 
the susceptibility is more pronounced towards the CEP. 
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T [MeV] T [MeV] 

FIG. 8: Entropy density normalized to the Stefan-Boltzmann value for a constant Tq — 208 MeV (left panel) and with 
/^-corrections (right panel). At the first-order phase transition the entropy density jumps. 




FIG. 9: Quark number density n q for two different chemical potentials as a function of the temperature for a constant To = 208 
MeV (left panel) and with /i-corrections Tb(/i) (right panel). The dashed lines are the corresponding Stefan-Boltzmann densities. 
The curve for fi = 290 MeV starts at the critical temperatures of the first-order phase transition. 



VI. SUMMARY AND CONCLUSIONS 

In the present work we have studied the Polyakov- 
extended quark-meson model beyond mean-field approx- 
imation. The quark-meson fluctuations to the matter 
sector are included within a functional renormalization 
group approach. In turn, the quark-meson fluctuations to 
the gluonic sector are estimated by a perturbative com- 
putation as suggested in Q . Interestingly, the latter pro- 
cedure has been recently confirmed by a full dynamical 
QCD computation [jj. fioj . 

The validity range of the present model includes the 
temperature regime about the crossover temperature 
T c w 200 MeV up to medium quark chemical potential of 
fi q w 100 — 200 MeV. At low temperature and large chem- 
ical potential the model suffers from missing inclusion of 
baryonic degrees of freedom and further resonances. 

However, the quantum fluctuations lead to strong mod- 
ifications of the phase structure. Such modifications have 
already been observed in FRG studies of the quark- meson 
model. One of the prominent effects is a shrinking of the 
size of the critical region around the critical endpoint in 
the phase diagram [421 ] . Furthermore, quantum fluctua- 



tions push the potential critical endpoint towards lower 
temperatures and larger chemical potential, see, e.g., f27j . 

Our computation entails a small likelihood for a critical 
point with /zr/T « 1 — 2 as predicted by some recent lat- 
tice studies 43 . Indeed, this is in accordance with the ar- 
guments of [44 | as well as with recent developments con- 
cerning the convergence of the Taylor expansion about 
vanishing chemical potential [45j]. The latter also com- 
plicates the extraction of the location of the critical point 
from the analysis of higher moments such as investigated 
in (46|. Thus, an extension of the present computations 
beyond mean field towards the low temperature high den- 
sity regime will provide valuable information. 

Further results concern thermodynamical quantities 
such as the pressure and the density. At vanishing den- 
sity these quantities agree well with related lattice pre- 
dictions and, in particular, show the correct behavior in 
the transition regime from mesonic degrees of freedom to 
the quark-gluon plasma regime. 

In a forthcoming publication [371 ] we will present a full 



solution including novel algorithmic differentiation tech- 
niques [36| for the flow Eq. (|2"Tj) which also allows us to 
discuss the systematics of the present approach towards 
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QCD, as well as studying further thermodynamic quan- 
tities. We also plan to extend the present study to the 
2 + 1 flavor well as tightening the relation to full 

dynamical QCD. 
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